Antiferromagnetic to superconducting phase transition in the hole- and 
electron-doped Hubbard model at zero temperature 
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The competition between d-wave superconductivity (SC) and antiferromagnetism (AF) in the 
high- Tc cuprates is investigated by studying the hole- and electron-doped two-dimensional Hubbard 
model with a recently proposed variational quantum-cluster theory. The approach is shown to 
provide a thermodynamically consistent determination of the particle number, provided that an 
overall shift of the on-site energies is treated as a variational parameter. The consequences for 
the single-particle excitation spectra and for the phase diagram are explored. By comparing the 
single-particle spectra with quantum Monte-Carlo (QMC) and experimental data, we verify that the 
low-energy excitations in a strongly-correlated electronic system are described appropriately. The 
cluster calculations also reproduce the overall ground-state phase diagram of the high-temperature 
superconductors. In particular, they include salient features such as the enhanced robustness of the 
antiferromagnetic state as a function of electron doping and the tendency towards phase separation 
into a mixed antiferromagnetic-superconducting phase at low-doping and a pure superconducting 
phase at high (both hole and electron) doping. 

PACS numbers: 71.10.-w, 74.20.-z, 75.10.-b 



I. INTRODUCTION 

The theory of high-temperature superconductivity re- 
mains one of the most challenging problems in solid- 
state physics. In many metallic compounds, over a 
broad range of compositions and temperatures, only two 
phases are encountered, i.e. the normal (Fermi liquid) 
and a magnetic or superconducting phase. In sharp 
contrast, the high-temperature superconductors (HTSC) 
and other strongly correlated electron systems, such as 
heavy-fermion and a variety of transition-metal oxide 
systems, exhibit many ordered phases, which appear 
to compete and sometimes coexisti In the HTSC, be- 
sides ordered antiferromagnetic (AF) and superconduct- 
ing (SC) phases, compelling evidence exists for charge- 
and spin- "stripe" phases and phases with coexisting SC 
and AF orderiS The common denominator and the un- 
derlying reason for these competing orders is certainly 
the presence of strong electronic correlations and Mott- 
Hubbard physics: The interplay between kinetic-energy 
and Coulomb-correlation effects induces an extreme sen- 
sitivity to external parameters (doping, temperature, 
pressure, etc.) and a rather difficult to predict "out- 
come" , i.e. the characteristic low-energy excitations and 
the phase diagram at low temperatures. The central chal- 
lenge in the field of high- Tc superconductivity is, there- 
fore, the connection of the (known) microscopic interac- 
tions at the level of electrons and ions, which are at high 
energy (^eV) and temperature T, with the "emerging 
phenomena" at T = 0, i.e. competing and nearly degen- 
erate orders. Ideally, one should employ a systematic 
renormalization-group approach to integrate out the ir- 
relevant degrees of freedom and, thereby, correctly bridge 
high to low energies and eventually go to T = 0. It is, 
however, by no means obvious how to do this when strong 



correlations are present, such as in the HTSC. 

In this context cluster techniques, which systematically 
approach the infinite-size (low-energy) limit, provide a 
powerful alternative^*^ In this paper, we will discuss 
and apply a variational-cluster approach (VCA) which 
was proposed recentlyi^ It is based on the self-energy- 
functional theory (SFT);^ which provides a general vari- 
ational scheme to use dynamical information from an ex- 
actly solvable "reference system" (in our case, an isolated 
cluster) to go to the infinite-size lattice fermion problem 
at low temperatures and at T = 0, in particular. In 
our earlier worki this scheme was formulated to study 
phases with spontaneously broken symmetry. For the 
cluster sizes used, it was shown that the VCA correctly 
reproduces long-range AF order for the two-dimensional 
{2D) Hubbard model and the absence of this order in one 
dimension. This non-trivial "test" implies that the VCA 
goes well beyond ordinary mean-field theory. 

Another crucial test is provided by the dynamical in- 
formation contained in the one-particle Green's function 
G. Compared to variational schemes based on wave 
functions an important advantage of the VCA consists 
in the fact that it quite naturally gives the one-electron 
Green's function G. For the 2D Hubbard model, it was 
recently demonstrated^ that the VCA, with the lattice 
tiled by (vTO x \/lO) clusters, correctly reproduces low- 
temperature quantum Monte-Carlo (QMC) data, in par- 
ticular, the coherent and incoherent "bands" experimen- 
tally known from ARPES dataiS 

These tests provide also the foundation for attacking 
the question whether the "minimal" microscopic model, 
namely the 2D one-band Hubbard model, reproduces the 
essential features of the electron- and hole-doped HTSC 
phase diagram. We will not go into a lengthy discussion 
of what interactions should be retained at the electron- 
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ion level. But, when choosing the 2D one-band Hubbard 
modelfii i.e. 

H = Y1 ^v4a(^j<y + UY1 "*T"a ' (1) 

ija i 

where Cjo-, c|g. are the usual annihilation and creation 
operators, tij denote the hopping-matrix elements, ni„ = 

cj^Cia the density at site i with spin cr =|, 4 and U the lo- 
cal Coulomb repulsion, one has introduced gross simplifi- 
cations, leaving out other orbital (e.g. p) degrees of free- 
dom, long-range Coulomb interaction, electron-phonon 
coupling, etc. Nevertheless, this choice appears to be le- 
gitimate, last not least in view of the amazing agreement 
achieved between numerical simulations and experimen- 
tal results for the normal-state properties of the cuprates 
(see, for example, R,efs. . 

The ground-state phase diagram of the model was re- 
cently investigated using the VCA by Senechal et al>i^ 
and by two of us^ There are important technical dif- 
ferences, but the "upshot" of the two works is as fol- 
lows: For the cluster sizes used in the VCA, the T = 
phase diagram of the Hubbard model (including hopping 
terms up to second or third-nearest neighbors) turns out 
to be qualitatively similar to that of the electron- and 
hole-doped cuprates. The model correctly describes the 
overall phase diagram, such as the occurrence of the AF 
and SC phases and predicts the corresponding doping 
ranges in qualitative agreement with the experiments for 
the cuprate materials. 

The present paper has several purposes: First, we 
would like to stress that for an application of the VCA to 
the high-Tc problem it is of crucial importance to treat 
the on-site energies in the reference system as variational 
parameters. We will show that this ensures a thermody- 
namically consistent determination of the average par- 
ticle number. Compared to the study of Ref. [ij this 
represents an important methodical extension. On the 
other hand, without the variational optimization of the 
on-site energy, one has to tolerate an inconsistency in 
the determination of the average particle number. The 
effects of this error shall be demonstrated by model cal- 
culations. The issue of thermodynamic consistency is 
also discussed for the off-diagonal elements of the one- 
particle density matrix and, in case of spontaneous sym- 
metry breaking, for the respective order parameter. It is 
interesting to note that there are no such problems in the 
dynamical mean- field theory (DMFT)iSiJI and its cluster 
extensions liSiiS Here the on-site energies are kept at their 
"physical" values from the very beginning rather than be- 
ing determined from the self-consistency condition. The 
(cellular) DMFT, however, can be considered as a spe- 
cial approximation within the general SFT frameworki&i 
Hence, the question arises why a fixed on-site energy does 
not spoil thermodynamic consistency in the case of (cel- 
lular) DMFT. It is interesting to note that, in case of 
interacting bosons, the issue of consistency has also been 
shown to be very important, recently^ 



Second, an accurate analysis of the behavior of the 
chemical potential as a function of the particle density 
close to the transition to a non-magnetic state, as well 
as a corresponding Maxwell construction, indicate the 
presence of an inhomogeneous ground state with macro- 
scopically large regions of low- and high-particle density. 
Using Lc — 4-site clusters, we can get a rough estimate 
of the ground-state phase diagram and investigate the 
instability of the homogeneous (AF, SC) phases against 
charge inhomogeneities. This represents an important 
complement to the work of Ref. [Ij, where larger cluster 
sizes up to Lc = 10 sites have been considered but with- 
out an appropriate analysis via a Maxwell construction. 

Finally, our numerical results give valuable insights 
into different questions of the high-Tc problem, as they 
provide direct access to the single-particle excitation 
spectrum in the strong-coupling regime at zero tempera- 
ture. This has to be contrasted with the work of Maier et 
aliSLS^ who have been able to treat clusters with Lc > 20 
sites within the dynamical cluster approximation (DCA) , 
but are restricted to finite temperatures, intermediate 
coupling U /t < A and imaginary time. A review compar- 
ing the application of different cluster (SFT, DCA, and 
cellular DMFT) as well as weak-coupling methods for the 
Hubbard model has appeared recentlyi^ 

Our paper is organized as follows: We start with a sum- 
mary of the central ideas of the SFT in general fSec. lII A|l . 
Thermodynamic consistency with respect to the average 
particle number is discussed in Sec. IIIBI Consistency 
with respect to the off-diagonal elements of the density 
matrix is addressed in Appendix ^ and the case of the 
C-DMFT is discussed in Appendix Our theoretical 
considerations are completed by describing some com- 
putational details of the VCA calculations in Sec. Ill CI 
Section IHII then presents the results for the T = spec- 
tral function A{k,u}) for electron- and hole-doping. We 
discuss how the characteristically different doping depen- 
dencies of the spectra give rise to different Fermi-surface 
evolutions upon doping. These Fermi-surface evolutions 
can then be tied up with the characteristic differences 
in the electron- and hole-doped phase diagrams, such as 
the enhanced robustness of the AF order in the electron- 
doped case. The ground-state phase diagram will be pre- 
sented and discussed in detail in Sec. IIVI Finally, Sec.lVI 
contains our main conclusions and a summary. 



II. THEORETICAL BACKGROUND 

A. Self-energy-functional theory 

The central idea of the self-energy-functional theory 
(SFT)^ is to make use of the universality of the Luttinger- 
Ward functional <^u[G^ or of its Legendre transform 
Fj7[S]: For a system with Hamiltonian H = Ho^t) + 
Hi{U), where t are the one-particle and U the interac- 
tion parameters, the functional dependence $i7[---] or 
-?[/[•■■] is independent of t. This universality is obvi- 
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ous as the Luttinger-Ward functional is defined via a 
skeleton-diagram expansion involving dressed propaga- 
tors and vertices onlyiSi 

Concentrating on the self-energy S instead of the 
single-particle Green's function G, the grand potential 
of the system at temperature T and chemical potential 
/i can be written functional of S; 



r!t.[/[S] =Trln(Go 



-1 
t 



(2) 



where Gq,* = + /i — t) ^ is the free Green's function 
and Tr = r V , e'"^"°^tr with the usual trace tr and the 
Matsubara frequencies a;„ = (2?7 + l)7rr for integer n. At 
the physical self-energy S = St f/i the grand potential is 
stationary: 5Q.t,u\^t.u] ~ 0. 

Why is it more advantageous to express as a func- 
tional of the self-energy S rather than G? This has 
to do with the "short-range" character of S as a func- 
tion of its real-space coordinateSfSi which in general is 
due to the fact that S is qualitatively related to a dy- 
namically screened particle-particle interaction. For the 
Hubbard model, in particular, one believes that impor- 
tant effects are sufficiently accounted for by a localS^ or 
short-ranged2& self-energy, at least for high lattice dimen- 
sions. A local or short-ranged self-energy, however, can 
be well generated by a cluster of finite size, and for the 
subsequent optimization of the cluster trial self-energy, 
the self-energy functional ^ can be used. This con- 
cept allows for the construction of a class of conceptually 
clear and thermodynamically consistent approximations, 
including the dynamical mean-field theory (DMFT}iS*ii 
and a cluster extension of the DMFTii^ (see Ref. |23 for 
a detailed discussion). 

Due to the universality of Fu [S] , we have 



^!t^c/[S]=Ff/[S]+Trln(Go 



-1 
t' 



(3) 



for the self-energy functional of a so-called "reference sys- 
tem" , which is given by a Hamiltonian with the same in- 
teraction part U but modified one-particle parameters t': 
H' = Ho{t') + Hi{U). Although it has different micro- 
scopic parameters, the reference system is assumed to be 
in the same macroscopic state as the original system, so 
it has the same temperature T and the same chemical po- 
tential fj.. By a proper choice of its one-particle part, the 
problem posed by the reference system H' can be much 
simpler than the original problem posed by H, such that 
the self-energy of the reference system, St'.u, can be 
computed exactly within a certain subspace of parame- 
ters t' . Combining Eqs. (0) and Q, we can eliminate 
the functional Fu [S] . Inserting as a trial self-energy the 
self-energy of the reference system, then yields: 

nt.u['^t',u] = nt'.u+Tiln{Gal-'St'.uy^-Ti-lnGt',u, 

(4) 

where flt\u and Gt'jj are the grand potential and the 
Green's function of the reference system. Stationary 
points are obtained, and this is the approximation, by 



restricting the variation to the subspace of trial self- 
energies St'.t/: 



da 



t,ul-^t\u\ 
'dt' 







for t' = t 



(5) 



Varying the trial self-energy means to vary the one- 
particle parameters t' of the reference system. For fur- 
ther details of the approach see Ref. |^. 

Here, we will focus on the Hubbard model, Eq. 
as the original model given by H. Different possi- 
ble choices for H' and the corresponding systematics of 
approximations generated in this way are discussed in 
Ref. I23. In the following, we will concentrate on two 
cluster approaches: (i) the variational cluster approach 
(VCA)&i and (ii) the cellular DMFT (C-DMFT),i& The 
VGA can be seen as a variational generalization of the 
cluster-perturbation theoryiSSiSS It is obtained by parti- 
tioning the infinite lattice into disconnected (identical) 
clusters of Lc sites each and choosing H' to consist of 
the intra-cluster parts only, i.e. the inter-cluster hopping 
is switched off in H'. The C-DMFT is a cluster variant 
of the DMFT. In the context of the SFT, it is obtained in 
the same way as the VGA but with an additional coupling 
of each of the Lc correlated cluster sites to a continuous 
bath, i.e. to an infinite number of uncorrelated additional 
bath sites. The on-site energies of the bath sites as well 
as their coupling to the original sites are treated as vari- 
ational parameters. 



B. Consistent determination of the particle density 

Once a reference system is specified, one should, in 
principle, vary all one-particle parameters of H'. This 
procedure would give the optimal result but requires a 
search for a stationary point in a high-dimensional pa- 
rameter space. From a pragmatic point of view it is thus 
advisable to concentrate on a few parameters only which 
have to be selected by physical arguments. Here, we ar- 
gue that the variation of the on-site energies is impor- 
tant to achieve thermodynamic consistency with respect 
to the average particle number. In case of the Hubbard 
model, this means to consider the site-independent en- 
ergy e' = i-j as (one of the) variational parameter (s). 

The average particle number (N) can be calculated in 
two different ways: on the "zero-particle level" by dif- 
ferentiation of the grand potential with respect to the 
chemical potential ^: 



(6) 



and on the "one-particle level" by frequency integration 
of the one-particle excitation spectrum: 

/oo 
f{u;)Aua{u;)doj . (7) 
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Here, Aiia(ijj) = — ImGiio-('^ + *0+)/7r is the local (possi- 
bly spin-dependent) spectral density with G = Gt^jj for 
short, and /(w) = {exp{uj/T) + 1)^^ is the Fermi func- 
tion. For simplicity, we exclude non-coUinear magnetic 
states and assume all expectation values to be diagonal 
in the spin index. 

Thermodynamic consistency means that both ways 
of calculating {N) yield the same result. Since = 
VLt^u\^t' ,u] is the approximate SFT grand potential at 
the stationary point t' = t'^, and since the spectral 
density or, equivalently, the one-particle Green's func- 
tion is the approximate Green's function given by G = 
1/(Gq J — St^^[/)~^, the equivalence of ^ and O is by 
no means understood a priori. 

To prove thermodynamic consistency, we start from 
Eq. According to Eq. there is a twofold ^ de- 
pendence of f2 = rit,u\^t'^.u]'- (i) an explicit fi depen- 
dence due to the chemical potential in the free Green's 
function of the original model, G'^^l = uj + fi — t, and 
(ii) an implicit fi dependence due to the fi dependence 
of the self-energy Ht'^^u, the Green's function Gt'^^u and 
the grand potential Clf jj of the reference system: 



an 



(8) 



Note that the implicit fi dependence is due to the chemi- 
cal potential of the reference system which, by construc- 
tion, is in the same macroscopic state as the original sys- 
tem (with the same temperature T and the same chem- 
ical potential fi) as well as due to the fi dependence of 
the stationary point tg itself. This is a subtlety which, 
however, can be ignored since 



on dt' _ 



(9) 



for t' = t'g because of stationarity condition (Actu- 
ally, only those elements of t'^ show up a dependence 
that are treated as variational parameters. According to 
the chain rule, however, the derivative of f2 has to be 
performed just with respect to those elements, with a 
vanishing result due to the stationarity condition.) 

The self-energy, the Green's function and the grand 
potential of the reference system are defined as grand- 
canonical averages. Hence, their dependence due to 
the grand-canonical Hamiltonian Ti' = H' — fiN is (apart 
from the sign) the same as their dependence on e': 
d/dnim. = —d/de', etc. Consequently, we have: 



{N) 



an 



an 



(10) 



The first derivative is readily calculated: 

an a 



9^cx. a^icK. 



tV e'-""^rln— ^ ^ 

"Z" Go ,(«^n) - St/ t/(«W„) 



-1 



iuJn + ll-t - St'_t/(«W„) 



— / e"°^/(cc;) tr ^— — 

2-1:1 J Q Lu + - t - IJt'^^ui^^) 



duj . 



(11) 



Here, the contour C encloses the first-order poles of the 
Fermi function at ujn = (2n + 1)ttT in counter-clockwise 
direction. Using Cauchy's theorem, we can proceed to an 
integration over real frequencies. Inserting into Eq. l(TU|l . 
we get: 



1 /" 

(N) = Im / 

J-, 



f{uj) tr 



1 



t',i/ 



duj + 7— 7 



an 

(12) 
for t' = 

The first term on the right-hand side is just the expres- 
sion for the average particle number given by IJJ). The 
second term on the right-hand side vanishes provided that 
the variational condition JSJ is satisfied, i.e. provided that 
e' is included in the set of variational parameters. In this 
case one has thermodynamic consistency. If e' was not 
treated as a variational parameter but kept at the value 
given by the original system, e' = e, one would have a 
finite an/ae' in H12|) . and the two expressions l|njl and 
([7|) for the average particle number would yield different 
results. This completes the proof. 

Eq. ^ for the average particle number of the lattice 
model can be compared with 



{N)' 



an' _ ^ 



f{co)A^,Aco)dco, (13) 



which gives the average particle number of the reference 
system. Again, there are two ways to get (N)': either as 
the derivative of the reference system's grand potential 
n' = nt'^u or by frequency integration of the reference 
system's spectral density A'--^{uj) = (— l/7r)ImG-j^(a' -I- 
lO^). As the reference system is solved exactly, both 
ways must yield the same result. Note, however, that 
(N) ^ (TV)' in general. 

The above reasoning can straightforwardly be gener- 
alized to the off-diagonal (j j) elements of the one- 
particle density matrix {c[„Cja). This is discussed in Ap- 
pendix 

The effect of thermodynamic (in)consistency is illus- 
trated for the single-band Hubbard model of Eq. (pj in 
Fig- n The figure shows n = {ni„) as a function of ^ as 
obtained from VGA calculations described in Sec. Ill Cl be- 
low. Solid lines display the result obtained by frequency 
integration of the spectral density, Eq. {Tj); while dashed 
lines show the result of the numerical /i derivative of 17, 
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FIG. 1: Filling n = (riic,) as a function of chemical poten- 
tial ^ obtained by integration of the spectral density (Eq. Q , 
solid lines) and via the derivative of Q, (Eq. dashed lines). 
Results are obtained, respectively, by considering e' as a vari- 
ational parameter (a), and by setting e' = (b). Calculations 
for the electron-doped case with (7 = 8, nearest-neighbor hop- 
ping tnn = — 1 and next-nearest-neighbor hopping tnnn = 0.3. 



Eq. © . We compare the results obtained by considering 
e' as a variational parameter (a) , with the ones obtained 
by setting e' = (b). In case (a) the two curves coincide 
within numerical accuracy, as expected, while in case (b) 
a considerable discrepancy is observed. This discrepancy 
increases upon approaching the transition to the non- 
magnetic state, i.e. precisely in the interesting region, 
where it reaches about 5% doping. 

While thermodynamic consistency with respect to the 
particle number is an issue for most approximations 
within the SFT, there is one exception: the cellular 
DMFT. This is discussed in Appendix |b1 



C. Variational cluster approach (VGA) 

In detail, the calculation proceeds as follows: We con- 
sider the 2D single-band Hubbard model of Eq. ^ with 
nearest- and next-nearest {tnnn) neighbor hoppings 
and the Hubbard U . For the purpose of our qualitative 
description of the HTSC, it is sufficient to take parame- 
ters which are typical for both hole- and electron-doped 
high-Tc cuprates, i.e. tnnn/tnn = —0.3 and U/tnn = 8. 
The energy scale is set by choosing <„„ = —1. Differ- 
ent values of the parameters (e.g. a change of U or t„„„ 
within 30%) or a third-nearest neighbor hopping have 
been also incorporated, checked and found not to quali- 
tatively change our conclusions (for example, concerning 



the weak phase separation detected in the electron-doped 
case). 

The Hamiltonian of the reference system H' is given by 
a set of decoupled clusters of finite size. For an individual 
cluster, the Hamiltonian reads: 



n' — TJ' 

-f^cluster — ^Hub. 



(14) 



It consists of the Hubbard Hamiltonian confined to the 
finite cluster H^^y^ plus two symmetry-breaking terms 
(Weiss fields) iJ^F ^^'^ ^sc with 

^;F-/ikFE("''T-"a)e^'^''- (15) 



and 



^sc = ^sc ■^"('^^TCji + h.c.) 



(16) 



where h'^p is the strength of the staggered and h'^^ 
the strength of the nearest-neighbor c?-wave pairing field. 
Q ~ (tt, tt) is the AF wave vector, and rjij denotes the 
d-wave form factor which is non-vanishing for nearest- 
neighbor lattice sites only and is equal to +1 (—1) for 
Ri — Rj in X (y) direction. The sum in Eq. (|16|) is 
restricted to sites i and j belonging to the same clus- 
ter. According to the discussion in Sec. Ill Bl the site- 
independent energy e' ~ t'^^ shall be treated as a varia- 
tional parameter. It shows up in the local term 



(17) 



which is already included in H^^^^ . The optimization of 
e' has to be done simultaneously with the optimization 
of the parameters /i^p and hg^. 

Due to the optimization of the Weiss field strengths 
/i^p and h'gQ, one can account for spontaneous AF and 
d-wave SC symmetry breaking. The respective AF and 
SC order parameters, m and A, are defined as 



on 

(9/lAF 



A 



dh 



SC 



(18) 



in the limit hAF,hsc where hfi^p^hgc are the 
strengths of external physical staggered and pairing fields. 
These physical fields should not be confused with the fic- 
ticious Weiss fields with strengths /i^p and /ig^. Adding 
the respective physical field terms to the Hamiltonian H 
and performing the derivative with respect to /iAF,/isc 
of the SFT grand potential (at the respective optimal 
ficticious field strengths /iAF'^Sc)' yields m and A con- 
sistently with their representations (on the "one-particle 
level") as frequency integrals over the usual and anoma- 
lous one-particle spectral density. This consistency is 
shown in the same way as in Sec. Ill Bl for the average 
particle number and is a consequence of treating the fic- 
ticious fields as variational parameters. 

The quality of the approximation is decisively influ- 
enced by the cluster size used. On the one hand, for an 
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appropriate characterization of the phase transition (as 
considered in Sec. IV, below), one needs a sufficient ac- 
curacy in the grand potential ft. This accuracy is, first of 
all, determined by the requirement that the clusters cho- 
sen must be large enough to fully account for the "short- 
range" spatial dependence of the self-energy as has been 
discussed below Eq. |j2Jl. On the other hand, for a given 
cluster size, an as accurate as possible numerical evalua- 
tion has to be employed. 

With respect to the latter, we found it both conve- 
nient and accurate to evaluate the frequency integrals 
contained in the trace in Eq. Q not by numerical integra- 
tion (where the required accuracy is difficult to achieve) , 
but by converting these integrals to a sum over the poles 
of the Green's function (see Ref.l^for details). This, how- 
ever, requires the computation of all many-body eigen- 
states in a given sector of the cluster Hamiltonian H'^^^^^^^. 
(Eq. (|14() ') including symmetry-breaking fields. This tech- 
nically limits the cluster sizes to be considered. There- 
fore, in the present work, we have chosen an infinite lat- 
tice tiled with 2x2 clusters. 

Cluster consisting of 2 x 2 sites, as well as larger cluster 
sizes have recently been systematically studied by Kyung 
et al44 in their influence on various physical quantities 
of the 2D Hubbard model by means of the C-DMFT. 
Already the smallest, i.e. 2x2, cluster has been found to 
account for more than 95% of the "correlation effect" in 
the single-particle spectrum (for a precise definition see 
Eq. (16) of Ref . . This suggests that at least some of 
the relevant questions in a strongly correlated electron 
system, modelled by a Hubbard-type Hamiltonian, may 
be studied; even rather accurately, with a such a small 
reference cluster. 

Nevertheless, for a more accurate determination of the 
SC state and of the phase transition, larger clusters have 
to be considered eventually. Only in this way, phase- 
fluctuation effects between different d-wave Cooper pairs 
(which determine the spatial dependence of the self en- 
ergy in the SC state) can be accounted for. We expect, 
for example, that calculations for larger clusters would 
display a smaller SC gap than the one seen in Fig. ^ 
around (tt, 0) in the electron-doped case. On the other 
hand, we expect our results on the global phase diagram 
to be quite robust. This is corroborated by the only weak 
dependence on (moderate) model changes, such as the 
next-nearest hopping or the Hubbard interaction^i*i^ 



III. SINGLE-PARTICLE EXCITATIONS 

The evolution of the single-particle excitation spec- 
trum A{k,uj) as a function of hole doping, and, in par- 
ticular, the transition from a "small" Fermi surface (hole 
pockets around {tt/2, 7r/2)) to an LDA-like Fermi surface 
closed around (tt, tt) are key observations in unlocking the 
mystery of the cuprates. Earlier QMC calculations for 
the 2D Hubbard modeliiii^ found that the single-particle 
spectral weight A{k,uj) semi-quantitatively reproduces 



both the momentum (d^^-yi -symmetiy) and, in particu- 
lar, doping dependence of the "high-energy" pseudogap 
of the order of the exchange energy J ~ 200meV as found 
in photoemission experiments around (tt, 0). The corre- 
sponding QMC data are reproduced in Fig. [21 for com- 
parison with the results of the VCA shown in Fig. O for 
hole dopings x — O.OI and x = 0.05. 

We first discuss the QMC data. In the under-doped 
regime, the "pseudogap" feature near (tt, 0) moves to 
lower binding energy as doping is increased. At about 
X = 0.13 (full circle in Fig. |2J) the pseudogap vanishes in 
overall accordance with the experimental findings^ In 
the experiments, the "high-energy" pseudogap is identi- 
fied with the centroids of spectral weight near (tt, 0). 

However, the experimental ARPES spectra also dis- 
play a "low-energy" pseudogap in the normal (supercon- 
ducting) state above (below) Tc with energy ^ 20meV, 
inferred from the leading edge in the spectral density, 
which also opens up in the under-doped regime and van- 
ishes in the over-doped regime^ This empirical corre- 
lation between the disappearance of the "high-energy" 
pseudogap and the decrease of the SC gap and, therefore, 
pairing strength suggested already several years ago that 
the high-energy features at (tt, 0) (which are correctly 
described by the "high-temperature" QMC simulations) 
are closely related to the pairing interactioniiS 

It is clear, however, that there is a strong need to 
perform calculations at much lower temperatures, i.e. 
temperatures below Tc ~ 2QmeV. Only then the "low- 
energy" SC or normal state pseudogap can be detected 
and only then the question "where do holes enter first?" 
can be answered correctly. According to the results of the 
"high-temperature" (T = 0.33) QMC simulations shown 
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FIG. 2: Taken from Ref. H^ . The dispersion of the peaks in 
the single-particle spectral weight from (a) QMC simulations 
of the Hubbard model at the temperatures indicated and at 
dopings X = 0.01, x = 0.05 (peaks in A(k,uj) represented 
by error bars) and x — 0.13 (solid circle), (b) ARPES ex- 
periments from under-doped and optimally doped materials 
(peak centroids in A{k,uj)) after Ref. l30l 
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FIG. 3: Single-particle spectrum for the hole-doped {x = 0.01 
and X — 0.05) case at U — 8 = —1). For comparison 
with the QMC data, the next-nearest-neighbor hopping is set 



to t„ 
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FIG. 4: Single-particle spectrum for U — 8 and i„„„ = 0.3 
{tnn ~ —1) in the hole-doped (a) and in the electron-doped 
(b) case. Results are shown for dopings in the mixed AF-I-SC 
state (see Sec.|l3, i.e. for x = 0.015 (a) and x = 0.09 (b), 
respectively. 



in Fig. 121 for x = 0.01, holes enter first around (tTjTt). 
This can be understood by referring to the ideai^^ that 
a higher temperature T effectively acts as an increased 
doping which destroys the magnetic Brillouin zone. This 
allows holes to first enter into the "arc" of single-particle 
excitations spanned around Q ~ (tt, tt). 

However, a portion of the corresponding "large" Fermi 
surface seems to disappear already at the lower temper- 
ature T = 0.25. This is indicated in Fig. |2a) by the 
downturn of the quasi-particle-like band between (tt, 0) 
and (tt, tt). In the experiment (Fig. [2tb)), this behavior 
in the under-doped regime has been interpreted as the 
opening of a pscudogap in the underlying Fermi surface 
near the (tt, 0) to (tt, tt) Hnepi^ This opening is obviously 
intimately related to the above question, namely where 
doped holes first enter, i.e. to the possibility of hole pock- 
ets developing at very small temperatures and low dop- 
ings around (7r/2,7r/2) in the hole-doped case. 

In contrast to the finitc-T QMC result, the correspond- 
ing VGA calculations for T = (see Fig. ^ show that 
holes first go into the coherent band around (7r/2,7r/2) 
forming "hole pockets" consistent with experiments^^ 
Note that for the dopings considered in Fig. |31 the sys- 
tem is in a mixed AF-I-SC state with non-vanishing AF 
order parameter (see also Fig. Cj). Therefore, although 
the low-lying spectral weight is found around (7r/2, 7r/2), 
the spectrum still shows an AF gap at this wave vector. 



Related to this, it appears that even for x = 0.05 the 
quasi-particle-like dispersion does not cross the Fermi en- 
ergy along the nodal direction (between (0, 0) and (tt, tt)) 
in contrast to the experimental situation in Fig. ^jp. 
Although some (weak) accumulation of spectral weight 
close and above w = can be seen in Fig. we believe 
that this non-crossing is (partly) a shortcoming of our 
2x2 calculation and that the self-energy in the nodal 
direction requires larger clusters. 

To address this issue for more realistic model pa- 
rameters, i.e. including a ncxt-nearest-neighbor hopping 
tnrm = ~0.3i„,i, as wcU as to study the corresponding 
doping evolution in the electron-doped case, we discuss 
the VCA results for A{k,u!) shown in Fig.^J The calcu- 
lations have been performed for the hole- (a) and for the 
electron-doped system (b) in the mixed AF-I-SC phase 
(see Sec. II V|) . i.e. for x = 0.015 and x = 0.09, respec- 
tively. 

Let us concentrate on the very small (1.5%) hole 
doping first. In agreement with the corresponding 
experimentaifi in hole-doped cuprates, holes indeed first 
enter at (7r/2,7r/2). Although the system is in a mixed 
AF-fSC phase (see Sec. llVjl. the SC gap is zero at this 
nodal point, so that doping into nodal states, i.e. nor- 
mal metallic screening, apparently destabilizes the AF 
solution already for rather low doping values {x 0.03). 
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(0,51) 




FIG. 5: Evolution of the low-energy spectrum upon hole 
doping. Parameters are as in Fig. |1] The weight is obtained 
by integrating the low-energy j4(fc, tt;)-spectrum down to 0.2 
below the Fermi energy. 



This explains the phase diagram of Fig. {7\ 

A similar picture can be inferred by looking at the evo- 
lution of the Fermi surface as a function of doping. This 
can be extracted from Fig. [3 where the low-energy spec- 
tral weight is plotted in the Brillouin zone. In qualitative 
agreement with experiments, hole pockets start forming 
around (7r/2,7r/2) for low doping, while a large Fermi 
surface centered about (tt, tt) starts building up at higher 
doping. 

In contrast; in electron-doped systems, doped electrons 
initially form pocke ts around (7r,0) (see Figs.EJb) andEl 
sec also R.efs. ll4l3ll) . in agreement with experimentsi^SiS 
Fig. 1^ shows the spectral weight obtained by integrating 
A{k,uj) down to 0.2 below the Fermi energy (uj = 0). 
One has to be careful when interpreting Fig. since for 
X < 0.13 there is still an AF gap at the Fermi surface 
near (7r/2, 7r/2) besides the SC one near (tt, 0). There- 
fore, the excitation spectrum is completely gapped, and 
one has to go away from w = in order to find some 
weight. Since we integrate only in a small energy win- 
dow around the Fermi level, the scale in Fig.Elis an order 




10.025 



FIG. 6: The same as in Fig. |S| but for the electron-doped 
case. 



of magnitude smaller than in the hole-doped case, Fig. 
|31 for which the AF gap near (7r/2,7r/2) is shifted away 
from the Fermi surface. Nevertheless, one finds that the 
lowest-lying states in the electron-doped case are around 
(tt, 0). Here, the density of states is large and provides a 
large "reservoir" for electron doping. This, in combina- 
tion with the fact that the chemical potential lies in the 
SC gap (of the AF+SC phase), stabilizes the AF solution 
for a larger doping range than in the hole-doped case, al- 
lowing for the AF gap to decrease more gradually (see 
Sec. II V|) . Also in this case, the breakdown of the mag- 
netic solution occurs as soon as the chemical potential 
reaches the bottom of the band at {tt/2,tt/2). The ob- 
servation that doping into (7r/2,7r/2) generically makes 
the AF phase unstable, suggests that also the occurrence 
of phase separation, to be discussed in Sec. IIVI is rather 
generic and quite independent of model details. 

The low-energy spectral weight, as plotted for differ- 
ent dopings in Fig. El agrees quite well with the ARPES 
measurements of Armitage et alm^i^; doped electrons 
away from half- filling first form pockets around (tt, 0). 
At higher doping, Fermi-surface segments start develop- 
ing along the (tt, 0) — (0, tt) line. Finally, these segments 
connect and form a large Fermi surface around (tt, tt) near 
optimal doping. 



IV. GROUND-STATE PHASE DIAGRAM 

In this section the ground-state phase diagram and the 
AF to SC transition are presented and discussed (see also 
Refs. IT^ITsI) . In particular, we will focus on the similari- 
ties and differences between the hole- and electron-doped 
system as well as their relation with the evolution of the 
single-particle spectrum. 

In order to determine the T = phase diagram, we pro- 
ceed as described in Sec. Ill Cl i.e. two symmetry-breaking 
terms (Weiss fields) iJ^F ^^"^ -^sc ''^^'^ included and their 
respective strengths, ft-^p and hgQ are treated as vari- 
ational parameters in addition to e' . As discussed in 
Sec. Ill BI the use of the additional variational parameter 
e' is required in order to have a consistent determination 
of the particle density. 

The phase diagram for the Hubbard model with U = 8 
and next-nearest-neighbor hopping tnnn ~ 0.3 (we set 
tnn = ~1) as obtained from our calculations, is shown in 
Fig. 13 for the hole-doped and in Fig. [S] for the electron- 
doped case. In the upper part of each figure, the chemi- 
cal potential /i is plotted as a function of x. In the corre- 
sponding lower parts, we display the AF (m) and SC (A) 
order parameter as a function of doping x. Note that the 
order parameters m and A are different from the Weiss 
fields h'f^p and hgQ, respectively. Quite generally, how- 
ever, a nonvanishing stationary value for a Weiss field 
produces a nonvanishing order parameter, respectively, 
although the latter can be much smaller. 

Let us discuss hole doping first (see Fig. [Tj). For 
low dopings x we find a homogeneous symmetry-broken 



9 




0.05 0.1 0.15 0.2 



X 

FIG. 7: Antiferromagnetic and superconducting order pa- 
rameters, m and A, and chemical potential /x as functions of 
hole doping x. A and ^ are plotted for the AF+SC (green, 
Aaf+sc, MAf+sc) as well as for the pure SC homogeneous 
solutions (blue, Asc, Msc)- Note that A is scaled by a factor 
5 for convenience. For x < xi the system exhibits a coexis- 
tence of AF and d-wave SC order. Phase separation occurs 
between the doping levels xi and X2- For x > X2 pure d-wave 
SC is realized. In the phase separation region xi < x < X2, 
the homogeneous solution becomes unstable, and the system 
prefers to separate into a mixture of two densities correspond- 
ing to xi and X2. The chemical potential /ic is determined by 
the Maxwell construction shown in the upper figure. At fj,* 
the slope of the AF-I-SC solution changes sign. 



phase in which both the AF as well as the SC order pa- 
rameter m and A are non-zero. This corresponds to a 
phase where AF and SC order microscopically and co- 
herently coexist. A homogeneous phase with pure SC 
(to = and A > 0) is obtained for larger dopings. The 
behavior of m vs. x in Fig. [7| seems to suggest that the 
transition to the non-magnetic state is continuous (sec- 
ond order) as a function of doping. However, a glance 
at the non-monotonous behavior of the chemical poten- 
tial fj, plotted as a function of x in the upper part of the 
figure, indicates the occurrence of a charge instability. 
The system tends to separate into a hole-poor (xi) and 
a hole-rich (^2) phase. The two dopings xi and X2, as 
well as the chemical potential fic in the phase-separated 
region, are identified by the Maxwell construction shown 
in the upper part of Fig. [7| ^* is the point where the 
slope of fJ-{x) changes sign. 

Of course, this phase separation is obtained within a 
treatment that, in the present work, just considers for 
homogeneous order parameters only and that neglects 




X 



FIG. 8; Same as Fig. Q but for electron doping. Note the 
enhanced robustness of the AF state and the strongly reduced 
scale A/i = {fi* — fic) as compared to hole doping. 



surface effects. Furthermore, the VCA is a mean-field- 
type approach on length scales beyond the size of the 
individual cluster. Therefore, the above result has to be 
interpreted with care. We expect it to signal a tendency 
towards the formation of microscopic inhomogeneities, 
such as stripes, checkerboard patterns, etc^^t^^^ This 
tendency will be further studied in future work by consid- 
ering larger clusters and/or by allowing for a more general 
variational solution which explicitly describes stripe inho- 
mogeneities (for example, by considering coupled clusters 
with different dopings, as in Ref . . In this case, one 
might expect the phase transition to become more con- 
tinuous and phase separation to disappear eventually. 

Let us now discuss the electron-doped case and the 
similarities and differences as compared to hole-doping. 
The first observation is that, while the phase diagrams 
in Figs. 13 and |H1 arc qualitatively similar, the phase in 
which long-range AF order is realized is spreading to 
significantly larger doping values in the electron-doped 
case, in overall agreement with the experimental situa- 
tion. Fig. IHl shows that phase separation occurs in the 
electron-doped case as well, although the associated en- 
ergy scale A/i = (/i* — /ic) is smaller with respect to the 
hole-doped case by about an order of magnitude. In Ref. 
ITEI it is argued that this can be related to the different 
pseudogap and SC transition scales in hole- and electron- 
doped materials. This may give support to theories^^'^S 
which are based on the notion that fluctuations of com- 
peting phases, or of the related order parameters, are 
responsible for the pseudogap phenomenon. 
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FIG. 9: Q vs. h'p^p (with the two other parameters h'gQ and 
e' fixed at their stationary-point values) for difi'erent jj, in the 
hole-doped case. Parameters; U = 8, tnnn = 0.3 {tnn = — 1). 



Depending on the value of /i, there may be two solu- 
tions, an AF+SC and a pure SC one, corresponding to 
two stationary points of fi. In order to gain insight into 
the first-order transition between these two solutions, it is 
instructive to observe the behavior of the grand potential 
57 as a function of the variational parameter /i^p associ- 
ated with this transition for different ji. For simplicity, 
we consider here the hole-doped case only. Results are 
plotted in Fig.^ Note that /i^p is varied, while the other 
two variational parameters, h'^Q and e', are fixed at their 
values at the respective stationary point. At low doping, 
there is a single minimum of the grand potential at 
a finite value /i^p only (Fig. O /i = 1.35). In contrast, 
the paramagnetic (/i^F ~ 0) solution is given by a local 
maximum. Upon further doping, U, additionally develops 
a local minimum at vanishing AF variational parameter 
/i^p = (Fig. El M 1-30)- For fi < fic, the mini- 
mum at ft-AF — becomes lower than the one at finite 
/i^p indicating a first-order phase transition to a non- 
magnetic state (Fig. 1^ /i = 1.24). Eventually, the local 
maximum lying between the two minima merges with 
the minimum at finite /i^p s-nd the AF solution disap- 
pears (Fig. El n = 1.19). This disappearance just occurs 
when the chemical potential /.i enters the quasi-particle 
band around (7r/2,7r/2) {fi* = 1.22). 

Although Fig.Eldisplays the results for the hole-doped 
case, this behavior occurs qualitatively in the electron- 
doped case as well, although with a much smaller char- 
acteristic energy scale as discussed above. Again, these 
results reflect a qualitatively similar behavior for electron 
and hole doping. 



V. SUMMARY AND CONCLUSIONS 

By means of a recently developed quantum-cluster ap- 
proach, we have carried out a detailed analysis of the 
phase transition from the antiferromagnetic to the super- 
conducting phase in the Hubbard model at zero temper- 
ature. The main results concerning the nature of the AF 
to SC transition are summarized in Figs. [7| and |H1 At low 
dopings the AF phase actually mixes with a weak d-wave 
SC component. A similar coexistence phase is observed, 
for example, in PrCeCuO)2S and was also obtained in 
previous mean-field and cluster calculations>i^i2S4£ Upon 
further increasing doping, we find a transition to a pure 
d-wave SC phase. 

The phase-separation scenario, which is found for hole 
but also for electron doping within our approximation, 
should carefully be interpreted as a general tendency 
of the system to form microscopically inhomogeneous 
phases. Indeed, by allowing for a more general spatial de- 
pendence of the order parameter, the macroscopic phase 
separation will probably be replaced by other microscop- 
ically inhomogeneous phases, such as stripes, checker- 
board order, etc. In particular, this might be expected 
to be the case if long-range Coulomb interaction is taken 
into account additionallyi^^i^ii^ The situation is subtle 
in the electron-doped case. Here, in contrast to previous 
theoretical calculations, our results also suggest phase 
separation, although the corresponding energy scale (see 
Fig. [SJ is one order of magnitude smaller than in hole- 
doped compounds. 

Our results for the single-particle excitations presented 
in Sec. IIIII support the idea that it is, in fact, the single- 
particle spectrum which holds the key for understand- 
ing the qualitative differences seen in the ground-state 
phase diagram between the hole- and electron-doped 
cuprates. This concerns, for example, the robustness of 
the AF state in the electron-doped case. Here, electrons 
initially form pockets around (tt, 0) in accordance with 
experimentsi^ The fact that the density of states is large 
there, as well as the presence of a SC "gap" at this nodal 
point, stabilizes the AF state for a larger doping regime 
as compared to the hole-doped case. The breakdown 
of the magnetic solution appears as soon as "normal" 
metallic screening sets in, i.e. when fi touches the band 
at (7r/2,7r/2)i^ 

One should note that the variational cluster approach 
(VCA) is able to treat the fluctuations correctly up to 
the range of the cluster size only. Therefore, the question 
arises whether the SC solution we (and also othersi^i^Si^) 
obtain within the AF phase is a true long-range SC phase 
or whether it is only a signal of strong pairing fluctua- 
tions within the AF phase leading to a SC pseudogap. 
The latter hypothesis could be supported by the fact that 
results obtained with different cluster sizesi^ seem to in- 
dicate a size dependence of the SC order parameter, and 
by the fact that the SC order parameter is about a factor 
three smaller in the AF-I-SC phase than in the pure SC 
one. The presence or not of such a microscopic coexis- 
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tence phase may depend on material details. Certainly, 
our results suggest that the SC gap A (or pseudogap in 
the case of fluctuations) is important in order to stabilize 
the AF phase in electron-doped materials, as discussed 
in Sec. IIIll 



Similar VGA calculations have recently been carried 
out by Senechal et al>i^ using clusters up to 10 sites but 
without the variation of the on-site energies. These au- 
thors show that the single-band Hubbard model is suffi- 
cient to explain the different overall shapes of the phase 
diagrams for hole- and electron-doped cuprates. How- 
ever, their results seem to suggest that in the electron- 
doped case the AF to SC transition is continuous and 
associated with a quantum-critical point, in contrast to 
our results. This clearly shows the importance for a con- 
sistent determination of the average particle number. 

The recent substantial progress in relating the "high- 
energy" physics of the Hubbard model and its variants 
to the low-energy physics of competing phases is to a 
large extent due to the development of different quantum- 
cluster theories. Apart from the VGA, these are the 
cluster extensions of the DMFT, such as the cellular 
DMFT (C-DMFT) and the dynamical cluster approxi- 
mation (DC A). Kyung et ali^ have shown that some of 
the important problems in strongly correlated electron 
systems may be studied highly accurately using com- 
paratively small clusters. Maier et aliSLS^ performed a 
systematic cluster-size study of the 2D Hubbard model 
using rather large clusters (up to 26 sites). Converged 
results suggest a finite-T instability to d-wave SC state. 
Because of the QMC minus-sign problem, however, re- 
sults were limited to C/ = At where the typical correlation 
energy U and the magnetic energy scale of the HTSC is 
not yet achieved. On the other hand, the present VCA 
studies are clearly not yet converged with respect to the 
cluster size, as one can read off from, for example, the 
relatively large SC gap in the single-particle excitations, 
displayd in Fig. |5| An extension to larger clusters is, at 
least in principle, also possible within the VCA. This, 
however, necessarily implies the use of stochastic (QMC) 
methods as solvers for the cluster reference system. 



ONE-PARTICLE DENSITY 
MATRIX 



APPENDIX A: 



The reasoning in Sec. Ill Bl can straightforwardly be 
generalized to the off-diagonal (i ^ j) elements of the 
one-particle density matrix {cl^cja)- In this case, ther- 
modynamic consistency means that, for a selected pair 
of sites («, j), the derivative 

(Al) 



dU 



is equivalent with the integral 



f{uj) Ajm([j) duj 



(A2) 



where, for convenience, the hopping is formally assumed 
to be spin-dependent. 

Starting with Eq. (|Aip . we note that the tij„ depen- 
dence of f2 = ^t,u['^t' ,u] is due to the explicit tija- de- 
pendence of the free Green's function Go.t in Eq. 10}, and 
due to the implicit tij^ dependence of t'^: 



dt 



I] a 



TV e-"0+ ( ] 



dt' 



mi 

dt ■ 



(A3) 



Because of the stationarity of the grand potential, the 
second term on the r.h.s. can be ignored. We are thus left 
with the first term which, after transforming the Matsub- 
ara sum into an integration over real frequencies, exactly 
yields the desired expression (|A2|) . 

In case that not all elements of t' are treated as varia- 
tional parameters, the derivative of in the second term 
on the r.h.s. is a derivative with respect to those elements 
only. If, for the given pair of sites i and j, t[j^ is selected 
as a variational parameter, the contribution due to the 
implicit tyo- dependence vanishes again. If, on the other 
hand, t^^^ is not treated as a variational parameter but 
tij„ from the beginning, the second term 



fixed at t[^„ 



will give a finite contribution which spoils the thermody- 
namic consistency. 
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APPENDIX B: CELLULAR DMFT 

The C-DMFT is obtained by choosing as a refer- 
ence system disconnected clusters with a continuous non- 
interacting bath coupled to each cluster site. Carrying 
out the t' partial derivatives in Eq. and using Eq. 
for the SFT grand potential, we get the Euler equation 
in the form: 




= 0. (Bl) 
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Here i and j run over the sites of the original lattice, 
excluding bath sites. The one-particle bath parameters, 
namely the on-site energies of the bath sites and the hy- 
bridization of the bath sites with the correlated (physical) 
sites, have to be treated as (a continuous set of) varia- 
tional parameters. In the C-DMFT it is assumed that 
bath parameters can be found such that the first factor 
in ljBl|l vanishes: 



1 



r \ 



t'.U 



(B2) 



for arbitrary uj and for sites i and j belonging to the 
same cluster. Note that, by construction of the reference 
system, I],;jcr(^) and also dYiija{Lj) / dt' vanish if i and j 
belong to different clusters. Consequently, if bath param- 
eters can be found such that Eq. (|B2|) holds, the Euler 
equation ljBl|) will be satisfied, too: The self-energy func- 
tional is stationary at the C-DMFT self-energy. Eq. HB2|) 
is jus t the self-consistency equation of the C-DMFT (see 
Refs. lfillSl) . 

It is easy to see that {N) = {N)' within the C-DMFT: 
This simply follows by comparing Eqs. p2(l and for 
(iV) and {N)' and by using the C-DMFT self-consistency 
equation IIB2|) . 

Consider now the high-frequency expansions of the 
Green's function and of the self-energy of the reference 
system. Using G' = Gf u^^ = ^t'^c/. and Go = Go.t 
for short, we have (see Ref. |45|for high-frequency expan- 
sions in the Hubbard model): 



Here, {rii-a)' is the average occupation in the reference 
system. Using G'^\^^{ui) = [lo + iJ.)Sij — tij, inserting the 
expansions into the C-DMFT self-consistency equation 
(|B2|) . and expanding in powers of u;~^ once more, we 
immediately find: 



(B4) 



(B3) 



for (correlated) sites i and j within the same cluster, 
and, in particular, e' ~ e where e' = and s ~ ta. 
This means that, within the C-DMFT, consistency with 
respect to the particle number is assured by setting the 
on-site energies within a cluster of the reference system 
to their "physical" values. 



It goes without saying that a consistent determination 
of (TV) requires an exact solution of the reference system 
which, in case of the C-DMFT, is by no means trivial. 
An obvious idea is to replace the continuous bath by a 
few uncorrelated sites only to allow for an application of 
the Lanczos method (see Ref. for an example). This, 
however, immediately implies that the self-consistency 
equation (|B2p cannot be fulfilled exactly any longer and, 
strictly speaking, the determination of the particle num- 
ber becomes inconsistent unless s' is treated as a varia- 
tional parameter within the SET framework. 
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